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Abstract 

We  have  made  significant  progress  both  in  design  and  fabrication  of  Si-based  laser 
device  structures.  We  have  demonstrated  several  key  parts  of  the  laser  structure  including 
a  SiGe  waveguide  with  good  optical  and  electrical  properties  for  both  light  and  carriers 
confinement.  SiGe  resonant  cavity  with  DBRs  to  serve  as  Si-based  laser  cavity  was 
design  and  fabricated,  which  is,  to  the  best  of  our  knowledge,  the  first  SiGe  based 
resonant  cavity  with  DBRs  reported  thus  far.  In  addition,  we  have  carried  out  theoretical 
calculation  for  coupled  SAQD  systems  and  formulated  the  distinct  advantage  of  a 
coupled  SAQD  system  for  laser  applications.  In  conjunction  with  ARO  program 
DAAD19-01-1-0532,  we  have  achieved  good  control  and  the  mapping  of  of  the 
parameter  space  for  fabricating  InAs  quantum  dots  on  Si  substrate.  From  both 
experiments  and  theoretical  calculation,  we  have  proved  the  Volmer- Weber  growth  mode 
for  InAs  on  Si.  We  have  also  carried  out  detailed  transmission  electron  microscopy  (TEM) 
study  of  individual  dots  and  experimentally  identified  the  critical  dot  size  for  dislocation 
under  different  growth  conditions.  We  are  able  to  control  the  density  and  dot  size 
uniformity  of  InAs  SAQDs:  density  from  109  cm'2  to  1011  cm'2,  dot  diameter  from  9  nm 
to  30  nm,  dot  height  from  4nm  to  15nm  and  size  distribution  within  30%.  The  very  small 
critical  size  of  InAs  SAQD  for  dislocation  and  the  associated  extreme  quantum 
confinement  show  that  there  might  not  be  a  window  in  the  dot  size  for  which  electron- 
hole  pairs  can  be  confined  and  the  dots  are  coherently  strained. 
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Highlight  of  Research  Accomplishment 


1.  Laser  cavity  fabrication  and  measurement: 

1).  SiGe  waveguide  design  and  fabrication 

To  achieve  the  light  confinement  in  laser  device,  a  waveguide  structure  is 
needed  in  the  Si-based  quantum-dot  (QD)  laser.  SiGe  waveguide  is  a  good  choice  due 
to  its  low  loss,  good  coupling  with  fiber  and  flexibility  of  doping.  For  good 
performance,  the  parameters  of  the  waveguide  structure  need  to  be  optimized. 

Effective  refraction  index  method  and  mode  counting  technique  were  used  to 
get  the  phase  diagram  of  guided  modes.  In  these  diagrams,  a  safe  region  can  be  found 
in  which  only  single  TE  mode  and  TM  mode  are  supported  by  the  structure.  So  we 
can  choose  any  parameters  in  this  region  to  meet  the  other  requirements,  such  as 
critical  thickness  for  strained  layer  without  dislocations,  cross  section  matching  with 
single  mode  fiber.  Also,  field  distribution  in  the  waveguide  was  calculated  by  using 
preconditioned  block-iterative  eigensolvers  in  a  planewave  basis.  With  the  field 
distribution  (Fig.l),  optimized  position  for  active  layer  of  QDs  in  the  waveguide  can 
be  decided  at  the  highest  filed  intensity  area. 


-8  -6  -4  -2  0  2  4  6  8 

x  (um) 


Fig.l  Cross  section  filed  distribution  in  the  rib  waveguide. 

SiGe  waveguide  structure  with  optimized  parameters  has  been  fabricated  by 
Molecular  Beam  Epitaxy  (MBE).  Fig.  2  shows  the  propagation  loss  of  the  waveguide 
is  0.88  dB/cm,  which  should  be  compared  with  typical  waveguide  loss  in 
semiconductor  laser  structures  of  few  dB/cm.  It  should  be  pointed  out  that  the 
structure  we  measured  does  have  a  heavily  doped  capping  layer  necessary  for 
efficient  current  injection.  Fig.  3  is  the  I-V  characteristics  of  the  waveguide  structure. 
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To  understand  the  influence  of  metal  layer,  which  will  be  exist  as  electrode  in 
the  laser  device,  we  carried  out  a  study  of  waveguide  with  a  metal  layer  on  the  top.  A 
lpm  thick  silicon  spacer  (also  used  in  the  DBR  cavity)  is  used  to  introduce  better 
confinement  of  light  inside  the  SiGe  layer,  which  will  decrease  the  free  carrier 
absorption  in  the  metal  layer.  Also  such  a  layer  will  make  the  fully  strained  SiGe 
layer  more  stable.  The  propagation  loss  measure  (Fig.  4)  showed  such  a  structure  has 
a  loss  of  5.0dB/cm,  which  is  very  low  for  a  metal  cladding  waveguide  and  will  not 
cause  any  trouble  in  the  laser  cavity  applications. 


Length  (mm) 

Fig.4  Propagation  loss  measurements  of  SiGe  waveguide  with  300nm  A1 
layer.  The  waveguide  structure  are  shown  as  inserts. 


2).  Distributed  Bragg  reflector  (DBR)  laser  cavity 

A  high  quality  laser  cavity  for  the  lateral  mode  control  is  very  important  for 
the  single  mode  operation  of  laser  devices.  Transfer  matrix  method  with  an  additional 
gain  term  was  used  to  simulate  the  stimulated  light  emission  in  small  signal  condition. 
Fig.  5(a)  shows  the  emission  spectrum  from  standard  laser  cavity  with  one  layer  of 
InAs  QDs.  We  can  find,  due  to  the  inhomogeneously  broadened  gain  spectrum  caused 
by  size  fluctuation  of  the  quantum  dots,  there  is  a  large  number  of  longitudinal  modes 
even  in  a  very  short  cavity. 

To  suppress  unwished  longitudinal  modes,  the  distributed  Bragg  reflector 
(DBR)  laser  structure  was  used.  Phase  shift  condition  of  the  cavity  was  also 
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considered  in  the  design  to  obtain  more  suppression.  Fig.  5(b)  shows  the  stimulated 
emission  spectrum  from  such  a  structure.  Number  of  longitudinal  modes  is  greatly 
reduced,  so  single  mode  operation  can  be  easily  achieved. 


Mum) 


Fig. 5  Stimulated  emission  from  (a)  standard  cavity  of  l=14.4mm  (b)  DBR 
structure. 


Combining  the  grating  pattern  and  our  optimized  SiGe  rib  waveguide,  we 
fabricated  the  SiGe  resonant  cavity  for  the  laser  application.  Figure  6  is  a  SEM 
picture  of  the  resonant  cavity.  The  height  of  the  rib  is  0.7  pm,  the  etching  thickness  of 
the  grating  pattern  on  the  waveguide  is  0.3  pm,  the  length  of  the  cavity  is  21.84  pm, 
and  the  DBRs  at  both  side  of  the  cavity  include  1000  periods  grating,  where  each 
period  is  0.56  pm. 

To  characterize  the  resonant  cavity,  we  measure  the  transmission  spectrum  of 
the  resonant  cavity,  which  is  shown  in  Fig.  7.  Only  one  resonant  mode  centered  at 
1315.5  nm  can  be  found  in  the  transmission  spectrum,  which  shows  such  a  cavity  is 
very  suitable  for  single  mode  oscillation  of  SAQD  lasers.  From  the  transmission 
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spectrum,  we  can  also  calculate  the  cavity  quality  factor  0  =  317  and  photon  lifetime 
tc  =  0.22  ps,  which  can  be  further  improved  by  increasing  the  number  of  periods  in 
the  DBRs  and  the  etching  depth  of  the  trenches  in  the  DBR  pattern. 


8.0kV  1 2.5mm  x8.00k  SE(M)  3/23/04  1 5:01  5.00um 


Fig.  6  Top  view  SEM  photograph  of  the  SiGe  waveguide  with  DBRs. 


Fig.  7  Transmission  spectrum  of  the  SiGe  resonant  cavity  with  DBRs. 
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3).  Coupling  effect  on  the  laser  performance 


One  of  the  major  challenges  in  the  field  of  quantum  dot  laser  is  the  size 
ununiformity  of  quantum  dots,  which  will  cause  a  broad  gain  peak  and  therefore 
decrease  the  efficiency  of  laser  and  introduce  multi-mode  laser  operation.  A  lot  of 
effort  has  been  put  into  the  growth  of  highly  uniform  quantum  dots,  but  till  now  it  is 
still  very  difficult  to  get  a  very  uniform  quantum  dots.  In  this  report  we  will  show 
with  coupling  of  the  quantum  dots,  the  broadening  of  laser  spectrum  can  be  reduced 
and  the  mode  power  of  laser  can  be  increased.  And  a  system  with  poorer  size 
uniformity  will  benefit  more  from  the  coupling,  which  makes  coupling  even  more 
helpful  for  reducing  the  growth  difficulty. 

Fig.  8  shows  the  calculated  spectra  of  quantum  dot  laser  with  coupled  and 
uncoupled  quantum  dots.  The  calculation  is  done  by  solving  coupled  rate  equation  of 
such  a  quantum  dot  laser  structure.  In  the  results,  the  spectral  width  of  the  laser  output 
of  a  laser  structure  with  coupled  quantum  dots  is  significant  narrower  than  that  of  a 
system  with  uncoupled  quantum  dots.  This  phenomenon  can  be  understood  by 
considering  the  group  of  dots  that  resonant  with  the  most  intense  mode.  Because  of 
the  higher  stimulated  recombination  rate,  they  serve  as  a  fast  sink  for  excess  carriers. 
Coupling  allows  the  excess  carriers  from  the  dots  with  slower  recombination  rate  to 
migrate  into  the  sinks  leading  to  a  significantly  enhanced  emission  into  the  strongly 
resonant  mode  and  in  the  meantime  reducing  the  emission  into  the  other  modes.  The 
saturation  output  power  of  the  strongest  mode  of  the  laser  structure  with  different 
inhomogeneous  and  homogeneous  broadening  and  different  number  of  quantum  dots 
layers  are  also  calculated  using  coupled  rate  equations.  The  results  are  shown  in  Fig. 
9,  which  clearly  shows  that  the  enhancement  of  mode  power  by  coupling  effect 
increases  with  coupled  layers.  Another  important  thing  is  the  enhancement  due  to  the 
coupling  is  larger  in  the  system  with  poorer  dot  size  uniformity,  so  the  restriction  of 
dot  size  uniformity  can  be  loosed,  which  will  reduce  the  growth  difficulty. 


FIG.  8  Calculated  light  emission  spectra  of  quantum  dot  laser  with  coupled  and  uncoupled 

quantum  dots. 
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FIG.  9  Calculated  enhancement  of  single  mode  output  power  of  quantum  dot  laser  with  FWHM 
of  inhomogeneous  broadening  20meV  (circles)  and  10meV(up  triangles). 
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2.  III-V  quantum  dots  on  Si  (001): 

1).  Fabrication  of  InAs  SAQDs  on  Si  (001):  Comparing  to  Ge  SAQDs  on  Si,  the 
considerations  involved  in  the  growth  of  InAs  SAQD  on  Si  is  quite  a  bit  more  complex. 
Because  there  is  only  one  group  in  the  world  headed  by  Prof.  Bimberg  and  Prof. 

Gosele  of  Germany  that  is  doing  serious  research  in  the  fabrication  of  InAs  SAQD  on 
Si,  there  isn’t  much  information  available  in  the  literature.  For  example,  the  most 
important  question  of  whether  the  growth  is  in  Stranski-Krastanov  mode  (a  wetting 
layer  followed  by  3D  growth)  or  Vomer-Weber  mode  (non-wetting  3D  growth)  is  not 
known  for  sure.  Furthermore,  the  effect  of  As  over-pressure  on  the  dot  nucleation  and 
shape  is  not  established,  although  we  have  been  growing  in  the  As  rich  regime.  Finally, 
the  critical  dot  size  for  dislocation  has  not  been  determined. 

For  InAs  coverage  dependence  study 

Plan-view  TEM  images  of  InAs  dots  with  different  InAs  nominal  coverage  are  shown 
in  Figure  10.  The  density  and  size  dependence  on  InAs  coverage  are  shown  in  Figure 
1 1 .  TEM  results  indicate  that  dot  formation  begins  at  very  early  stage  of  the  growth  (< 
0.3  ML).  At  0.3  ML  of  InAs  coverage,  the  dot  density  is  approximately  1.7xl010  cm'2 
indicating  Yolmer-Weber  (VW)  growth  mode. 


Figure  10.  Bright-field  plan-view  TEM  of  InAs  dots  with  various  InAs  nominal  coverage  (a)  0.3  ML,  (b) 
0.7  ML,  (c)  1 .0  ML,  (d)  1 .4  ML. 
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Figure  1 1.  (a)  Plot  of  InAs  dot  density  dependence  on  the  nominal  coverage  and  (b)  Plot  of  InAs  dot  size 
dependence  on  the  nominal  coverage. 


Figure  12  shows  the  cross-section  TEM  of  typical  InAs  dots  with  the  different 
coverage.  With  increasing  InAs  coverage,  not  only  dot  size  but  also  shape  of  the  dots 
change.  Most  dots  have  an  irregular  shape  at  0.3  ML  coverage  that  changes  to  highly 
regular  huts  with  (11 1)  facets  at  0.7  ML.  At  1.0  ML,  portion  of  dots  evolves  into  dome 
shape  with  no  clear  facets.  At  1.4  ML,  nearly  100%  of  the  dots  are  dome  shape. 


(a)  (b) 


Figure  12.  Bright-field  cross-section  TEM  of  typical  InAs  dots  with  different  InAs  nominal  coverage:  (a) 
0.3  ML,  (b)  0.7  ML,  (c)  1 .0  ML,  (d)  1.4  ML. 

For  arsenic  pressure  dependence  study 
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Density  (cm  ) 


The  strong  dependence  of  dot  density  on  As  BEP  is  in  clear  contrast  to  its  weak 
dependence  on  InAs  coverage,  as  shown  in  Figure  13(a).  InAs  dot  density  increases  in 
a  near  exponential  fashion  from  4.3xl09  to  1.8xlOn  cm'2  as  As  BEP  is  decreased  from 
9.2x10'  to  1.2x10'  torr.  InAs  dot  size  and  volume  change  with  decreasing  As  BEP,  as 
shown  in  Fig.  13(b).  The  total  dot  volume  decreases  at  the  low  As  BEP.  This  indicates 
that  there  are  2-dimensional  InAs  structures  formed  between  dots  at  low  As  BEP. 


Figure  13.  (a)  Plot  of  InAs  dot  density  dependence  on  arsenic  BEP,  (b)  Plot  of  InAs  dot  size  and  volume 
dependence  on  arsenic  BEP. 

The  well-reconstructed  2x1  surface  completely  disappears  after  InAs  growth  at  low  As 
BEP,  indicating  a  highly  “defective”  As  monolayer,  i.e.  with  high  density  of  As 
vacancies,  as  shown  in  Figure  14.  The  fact  that  high  InAs  dot  density  is  obtained  only 
at  low  As  BEP  supports  the  proposed  mechanism. 
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Figure  14.  RHEED  pattern  of  InAs  growth  with  different  arsenic  BEP:  a)  9.2xl0'6  torr,  b)1.4  xlO'6  Torr, 
c)  3.0xl0'7  Torr,  d)  1.2xl0'7  Torr. 


For  temperature  dependence  study  of  dot  formation, 

The  left  side  of  figure  15  shows  bright-field  plan- view  [001]  TEM  micrographs  of 
samples  grown  at  different  substrate  temperatures  increased  from  295  °C  for  a)  to  410 
°C  for  e).  The  right  side  of  figure  15  is  the  histogram  of  size  distribution  based  on  the 
plan-view  TEM.  Gaussian  curves  were  fitted  on  the  size  distribution  histogram  of 
figure  15.  It  is  clear  that  increasing  substrate  temperature  results  in  broadening  of  dot 
size  distribution.  Gaussian  peak  width,  referd  to  as  0.849  full-width  at  half  maximum 
(FWHM),  increases  from  8  nm  to  14  nm  when  substrate  temperature  increases  from 
295  °C  to  410  °C.  Coarsening  is  likely  to  be  responsible  for  the  broadening  of  dot  size 
distribution. 


Figure  16  shows  the  dependence  of  dot  density  on  substrate  temperature.  Dot  density 
is  nearly  independent  of  substrate  temperatures.  The  dot  density  at  the  onset  of 
nucleation  can  be  written  for  the  well-studied  quantum  dot  systems  such  as  Ge/Si  and 
InAs/GaAs  [\1 2]  as  follows[3]: 


n  oc 


/  F 


\Ds  J 


(1) 


where  F  is  the  growth  flux  coming  from  source,  Ds  is  the  surface  diffusion  coefficient 
which  increases  exponentially  (. Ds=Doexp(-Ea/ksT ))  with  increasing  temperature,  a  is  a 
constant  associated  mainly  to  the  critical  cluster  size,  a  increases  with  increasing 
cluster  size  and  saturates  at  unity  for  critical  cluster  size  larger  than  10  atoms.  [4]  For 
epitaxial  growth,  a  can  be  treated  as  a  constant.[14,15]  n  in  equation  (1)  is  the  initial  dot 
density.  This  simple  model  cannot  explain  the  experimental  observation  that  the  dot 
density  is  nearly  independent  of  temperature  within  the  range  of  295  °C  to  410  °C 
under  a  constant  arsenic  BEP.  This  suggests  that  the  nucleation  model  as  in  equation 
(1)  that  has  been  shown  to  be  valid  for  Ge/Si  and  InAs/GaAs  systems,  does  not  apply 
to  InAs  growth  on  Si  under  our  growth  conditions. 
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Figure  16.  Arrhenius  plot  of  InAs/Si  dot  density  (y  axis)  dependence  on  substrate  temperature  (x  axis). 

Optimization  of  InAs  quantum  dots  on  Si. 

High  quality,  narrow  size  distribution  and  high  density  dots  can  be  achieved  by  InAs 
growth  at  lower  temperature  and  lower  arsenic  BEP.  In  order  to  improve  crystal 
quality  and  avoid  the  indium  segregation  under  high  temperature  treatment,  the  dots 
growth  is  followed  by  annealing  at  higher  temperature  and  higher  arsenic  BEP.  The 
increase  of  RHEED  intensity  after  growth  and  annealing,  compared  to  samples 
grown  at  low  temperature  (below  250°C),  confirmed  the  improvement  of  crystalline 
quality.  Figure  17  (a)  shows  the  plan- view  TEM  image  of  a  sample  that  was  grown 
at  substrate  temperature  of  250°C  with  arsenic  and  indium  BEP  of  7.2  xlO'5  Pa  and 
7.6x1  O'6  Pa  and  annealed  at  substrate  temperature  of  410  °C  with  arsenic  BEP  of 
1 ,9x  10"4  Pa.  The  ramping  time  from  250  °C  to  410  °C  was  set  to  5  minutes.  After  the 
120  second  annealing  at  410°C,  the  sample  was  cooled  down  to  room  temperature  in 
approximately  10  minutes.  This  growth  procedure  resulted  in  QD  density  of  5.7xl010 
cm'2.  The  size  distribution  is  shown  in  figure  17  (b).  The  Gaussian  curve  fitting  of 
size  distribution  gives  a  narrow  peak  at  around  1 1  nm  with  Gaussian  peak  width  of  7 
nm.  The  improvement  in  size  distribution  and  the  high  density  are  attributed  to  low 
temperature  growth  and  low  arsenic  BEP,  respectively. 
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Figure  17.  InAs  dots  grown  on  Si  (100)  at  the  substrate  temperature  of  250°C  and  with  arsenic  BEP  of 
7.2xl0'5  torr  and  post-annealed  at  410°C  and  1.9xl0'4torr  arsenic  BEP:  a)  bright  field  plan-view  [001] 
TEM  micrograph  and  b)  Histogram  of  size  distribution  of  InAs  dots. 


InAs  dot  density  dependence  on  III-V  ratio: 

InAs  quantum  dots  were  grown  under  different  III/V  ratio.  The  dot  density  change 
with  III/V  ratio  is  shown  in  the  figure  18. 
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Figure  9.  dot  density  dependence  on  III/V  flux  ratio.  For  the  growth  rate  data  the  As  BEP  was  kept  at 
1.4xl0'6  Torr,  for  As  BEP  data  growth  rate  was  set  at  0.02  ML/s.  All  growth  was  done  at  320°C 


The  data  from  both  constant  As  BEP  and  In  BEP  show  the  same  trend  that  dot  density 
is  nearly  linearly  dependent  on  III/V  flux  ratio.  This  behavior  provides  a  way  to 
control  the  In  As  dot  density  in  a  very  large  range  from  109  cm'2  to  1011  cm'2. 

2).  Growth  Mode  of  III-V  on  IV  substrate 

From  the  growth  of  GaAs,  InAs  on  Si  and  GaAs  on  Ge,  we  conclude  that  they  are  all 
in  Volmer-Weber  growth  mode.  Fig.  19  showed  the  reflection  high-energy  electron 
diffraction  (RHEED)  for  the  early  stages  (below  1  ML)  of  GaAs  on  Si  (a),  InAs  on  Si 
(b),  GaAs  on  Ge  (c).  Fig.  19  (d)  and  (e)  showed  the  RHEED  pattern  of  Si  and  Ge 
surface  before  growth.  It  is  very  clear  that  the  3-dimensional  island  growth  at  very 
early  stage  (below  1  ML).  Therefore  all  of  the  growth  are  Volmer-Weber  growth 
without  2-dimensional  wetting  layer  growth.  Fig.  20  showed  the  AFM  pictures  of 
GaAs  0.7ML  on  Si  (a),  InAs  0.7  ML  on  Si  (b).  AFM  data  confirmed  3-D  Volmer- 
Weber  growth. 


17 


Fig.  19.  RHEED  images  for  different  growth,  (a)GaAs/Si  at  320°C,  (b)  InAs/Si  at  320°C,  (c)  InAs/Ge  at 
320°C,  (d)  Si  surface  before  growth,  (e)  Ge  surface  before  growth. 


(a) 


(b) 


Fig.  20.  1  xl  |_inr  AFM  of  GaAs  (a)  and  InAs  (b)  growth  on  Si  at  the  initial  stages  at  InAs/Si  at 

320°C  with  arsenic  BEP  1.4x  I  O'6  Torr. 


3).  First-principles  calculations  of  different  InAs  overlayers  on  Si(001) 

To  understand  the  mechanism  of  initial  stage  of  InAs  grown  on  silicon  (001)  surface, 
first-principle  total  energy  calculations  were  carried  out  with  the  plane-wave 
pseudopotential  method  within  the  generalized  gradient  approximation  (GGA).  A 
conjugate-gradient  algorithm  is  used  to  minimizing  the  total  energy  therefore  obtain 
the  atom  positions.  Fig.  21  shows  the  positions  of  atoms  in  different  systems  with 
different  InAs  overlayers  on  silicon  (001).  We  calculated  the  total  energies  of  systems 
with  different  thickness  of  InAs  on  silicon  substrate,  which  is  shown  in  Table  I. 
Comparing  the  total  energy  of  two  systems  with  the  same  number  of  atoms,  while  one 
system  has  half  surface  area  covered  with  single  As  layer  and  half  covered  with  As- 
In-As-In-As  layer  and  another  one  with  the  whole  surface  covered  with  As-In-As 
layer,  we  can  find  the  former  one  has  a  total  energy  of  -(63. 40+92. 36)/2  eY  =  -77.88 
eV,  which  is  lower  than  the  total  energy  (-77.82  eV)  of  system  with  a  uniform  As-In- 
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As  layer.  For  a  system  with  1/3  area  covered  with  As-In-As-In-As-In-As  and  2/3  with 
single  As,  the  total  energy  is  -(106.97/3+63.40x2/3)  eV  =  -77.92  eV,  which  is  even 
lower.  So  the  system  prefers  to  form  InAs  islands  with  As-finished  silicon  surface 
rather  than  uniform  InAs  layer  on  silicon  surface,  which  causes  a  3D  growth  mode  of 
InAs  on  silicon  even  in  the  very  initial  stage.  In  the  real  case,  since  the  island 
formation  in  3D  growth  will  partially  release  the  strain  energy  that  is  comes  from  the 
lattice  mismatch  between  InAs  and  Si,  the  3D  growth  will  be  even  more  energy 
favorable. 


-0.5  0.0  0  5  1  0  1.5  -0.5  0.0  0.5  1.0  1.5 

[1l0](aSi/sqrt(2))  [110]  (aat/sqrt<2)> 


Fig.  21.  Calculation  positions  of  atoms  in  (a)  Si(001)-As,  (b)  Si(001)-As-In-As  and  (c)  Si(001)-As-In- 
As-In-As. 
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TABLE  I.  Total  energies  of  systems  with  different  InAs  overlayers  on  silicon.  The  unit 
area  used  in  the  table  is  the  area  for  one  surface  dimmer. 


Structure 

Total  energy  (eV/unit  area) 

Si(001)-As 

-63.40 

Si(00 1 )- As-In- As 

-77.82 

Si(00 1  )-As-In- As-In- As 

-92.36 

Si(00 1 )- As-In- As-In- As-In- As 

-106.97 

4)  Critical  size  of  threading  dislocation-free  dots  at  different  growth  conditions: 

1.  Critical  size  dependence  on  arsenic  BEP 

Plastic  relaxation  in  strained  InAs  islands  grown  on  Si  (001)  by  MBE  as  a  function  of 
arsenic  beam  equivalent  pressure  is  investigated  using  transmission  electron 
microscopy.  Analyses  of  the  spacing  in  Moire  fringes  show  a  clear  trend  of  larger 
InAs  lattice  parameters  with  increasing  island  size  as  a  result  of  the  elastic  relaxation 
of  three-dimensional  structures.  Our  experimental  findings  show:  (1)  the  existence  of 
a  critical  lattice  parameter  (0.597  ±0.0 lnm)  instead  of  a  critical  island  size  that 
determines  the  onset  of  dislocation,  and  (2)  a  clear  dependence  of  the  critical  size  of 
InAs  clusters  on  arsenic  BEP  during  MBE  growths  from  1.2  x  10'  to  9.2  x  10'  Torr. 
We  postulate  that  these  results  are  the  consequences  of  the  As  BEP  dependent  surface 
energy  of  the  various  facets  and  the  associated  InAs  island  shape  change. 

Fig.  22  showed  the  dark-field  plan-view  TEM  images  of  InAs  dots  grown  under 
different  arsenic  BEP.  The  Moire  fringes  of  each  dot  can  be  accurately  determined 
from  images.  The  lattice  parameters  can  be  obtained  from  spacing  of  Moire 

c  /:  n 

fringes. [  ,  ]  The  dislocation  can  be  identified  from  the  terminated  Moire  fringe. [  ] 


Fig.  22  Dark  field  (g  =  220)  plan-view  TEM  micrographs  of  InAs  islands  on  Si  (001)  substrate  with  As 
pressure  equal  to  1.2xlO'7Torr  (A),  1.4xlO'5 6 7Torr  (B)  and  9.2xlO'6Torr  (C).  Samples  are  grown  at  320°C 

with  0.7ML  InAs  coverage 


5  P.D.  Miller,  C.P.  Liu,  W.L.  Henstrom,  J.M.  Gibson,  D.G.  Cahill,  Y.  Huang,  P.  Zhang,  T.I.  Kamins,  D.P. 
Basile,  R.S.  Williams,  Appl.  Phys.  Lett.  75  (1999)  46. 

6  P.  Chen,  Q.  Xie,  A.  Madhukar,  Li  Chen,  A.  Konkar,  J.  Vac.  Sci.  Technol.  B  12  (1994)  2568. 

7  “Electron  Microscopy  of  Thin  Crystals,”  P.  Hirsch,  A.  Howie,  R.B.  Nicholson,  D.W.Pashley,  and  M.J. 
Whelan,  (Robert  E.  Krieger,  New  York,  1977),  p.372. 
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The  size  distribution  of  both  dislocation-free  dots  and  dislocated  dots  was  shown  in  Fig. 
22.  It  is  very  clear  that  the  relative  percentage  of  relaxed  versus  strained  islands  increases 
monotonically  with  BEPas  (from  15%  of  dislocated  InAs  islands  for  the  lowest  As  BEP 
to  70%  for  the  highest  As  BEP) 
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Fig  23.  The  histograms  of  dislocation-free  (black)  and  dislocated  (gray)  islands  as  function  of  their 
base  sizes  corresponding  to  the  same  set  of  samples  as  in  Fig. 22. 
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The  critical  sizes  of  dislocation- free  dots  grown  under  different  arsenic  BEP  are  shown  in 
Fig.  24.  The  critical  base  size  dimensions  of  coherent  islands  increases  significantly 
(from  13  to  28  nm)  with  increasing  arsenic  BEP  from  1.2xl0"7  Torr  to  9.2x1 0"6  Torr. 
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Fig.  24  The  average  in-plane  lattice  parameter  of  InAs  islands  as  a  function  of  their  base  sizes  for  the  same 
set  of  samples  as  in  Fig. 22.  Dislocation  occurs  at  a  well  defined  a  =  0.597±0.01  nm  (as  indicated  by 
horizontal  lines)  and  is  relatively  independent  of  island  sizes. 
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2.  Critical  size  dependence  on  substrate  temperature 
Relaxation  of  InAs  islands  epitaxially  grown  by  molecular  beam  epitaxy  on  Si  (001)  as 
a  function  of  the  substrate  temperature  is  investigated  using  transmission  electron 
microscopy.  Analysis  of  the  spacing  in  Moire  fringes  shows  that  an  increase  in  the 
average  InAs  lattice  parameters  with  increasing  island  sizes  in  both  elastically  as  well 
as  plastically  relaxed  islands.  Using  this  method,  a  critical  InAs  islands  size  for 
dislocation  is  determined  to  be  13.5±0.1  nm,  with  a  corresponding  critical  InAs  lattice 
spacing  of  0.579±0.01nm  along  (Oil)  crystallographic  directions.  The  critical 
dimensions  are  independent  of  the  substrate  temperature.  Although  no  dislocations  can 
be  observed  in  islands  smaller  than  the  critical  size,  the  size  of  the  largest  dislocation- 
free  islands  does  increase  significantly  (from  ~17  nm  to  ~24  nm)  with  increasing 
substrate  temperature  (from  295°C  to  410°C). 


Fig.25  Dark  field  (g  =  220)  plan-view  TEM  micrographs  of  InAs  islands  on  Si  (001)  substrate  for  samples 
grown  with  As  pressure  equal  to  1.4xlO'6Torr  with  total  coverage  0.7ML  at  substrate  temperatures  295°C 
(A),  345°C  (B)  and  4 1 0°C  (C),  respectively. 

Fig.25  shows  typical  plan-view  TEM  micrographs  of  samples  grown  at  three  different 
temperatures:  295°C,  345°C  and  410°C,  respectively.  It  is  apparent  that  increase  in 
substrate  temperature  within  this  range  leads  to  a  slight  decrease  in  the  average  density  of 
InAs  islands,  owing  presumably  to  the  increased  surface  diffusion  of  InAs  and  the 
associated  coarsening.  There  is  indeed  a  clearly  broadening  of  InAs  island  size 
distribution  as  can  be  visible  by  comparing  Fig.  26  A,  B  and  C.  Calculated  values  of  «,s/ 
are  plotted  as  a  function  of  the  island  base  sizes  in  Fig.  27  A,  B  and  C.  Both  strained 
(closed  symbols)  and  relaxed  (open  symbols)  islands  are  considered.  The  ratio  of  relaxed 
to  strained  islands  increases  from  20  to  35%  with  an  increase  of  substrate  temperature 
from  295°C  to  410°C.  There  is  a  corresponding  increase  in  the  average  dots  size  as  can  be 
seen  from  Fig.26  A,  B  and  C.  An  important  feature  of  Fig.26  is  that  the  size  distribution 
of  coherent  islands  is  nearly  temperature  independent.  The  substrate  temperature  increase 
results  primarily  in  a  broadening  of  the  range  of  island  sizes  within  which  coherent  and 
dislocated  islands  coexist.  This  range  is  from  13.7  nm  -17.5  nm  from  samples  grown  at 
295°C.  It  increases  to  from  13.7  nm  -  24  nm  for  samples  grown  at  410°C. 
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Fig.  26  InAs  islands  size  distribution  that  corresponds  to  the  same  set  samples  as  in  Fig. 25,  showing 
number  of  coherent  (black)  and  incoherent  (gray)  islands  as  function  of  their  base  size. 
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Fig.  27  Plots  of  an  average  in  plane  lattice  parameter  of  InAs  islands  as  a  function  of  their  base  size  for  the 
same  set  samples  as  in  Fig.25.  Horizontal  line  corresponds  to  the  critical  island  lattice  parameter  that 
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indicates  the  upper  limit  for  the  coherent  islands.  Vertical  dashed  lines  show  regions  where  both  coherent 
and  incoherent  islands  exist. 


5)  InAs  growth  on  relaxed  SiGe  buffer  layer 

InAs  dots  were  grown  on  relaxed  SiGe  buffer  layer.  Fig.  28  shows  the  typical  sample 
structure.  Fig.  29  showed  AFM  images  of  Ge  dots  on  relaxed  SiGe  buffer  layer  and 
InAs  dots  on  relaxed  SiGe  buffer  layer.  It  is  very  interesting  to  notice  that  InAs  dots  do 
not  preferentially  nucleated  over  the  dislocation  network.  Comparing  to  Ge  dots  on 
relaxed  buffer  layers,  the  growth  mode  is  YW  instead  of  SK  and  the  misfit  strain  is 
increased  from  4%  to  11%  for  InAs  on  Si.  Which  one  of  these  two  factors  prevents  the 
InAs  dots  from  registering  to  the  underlying  dislocation  network  is  an  important 
question  to  answer.  We  plan  on  conducting  studies  using  GaAs  on  SiGe  to  single  out 
the  chemical  incompatibility  from  the  misfit  strain  in  an  effort  to  hunt  for  the  answer  to 
this  question. 


Si  capping  100  A 


Relaxed  Sio.s  Geo.2  800  A 


Fig. 28,  The  substrate  structure  with  dislocation  network. 


(a)  Ge/  Si  (b)  InAs/  Si  (c)  InAs/Si 


Fig.  29.  (a)  5x5prrf  AFM  image  of  Ge  quantum  dots  on  Si,  (b)  5x5pm2  AFM  image  of  InAs  quantum  dots 
on  Si,  (c)  3x3  pm2  AFM  image  of  InAs  quantum  dots  on  Si. 


6)  Contrast  between  InAs  and  GaAs  on  Si  (001) 

Both  InAs  and  GaAs  are  grown  on  Si  (001)  with  the  similar  growth  conditions.  A  clear 
difference  between  the  InAs  and  GaAs  growth  behavior  is  the  dot  density  dependence  on 
arsenic  BEP.  InAs  dot  density  strongly  depends  on  the  arsenic  BEP,  while  GaAs  dot 
density  is  nearly  independent  on  arsenic  BEP,  as  figure  30  shown.  This  suggests  that 
there  is  fundamental  physical  difference  between  the  InAs  and  GaAs  grown  on  Si  (001). 
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Figure  30.  The  blue  line  is  the  GaAs  dot  density  and  black  line  is  the  InAs  dot  density.  The  growth  rate  is 
0.01  ML/s  for  GaAs  and  0.02  ML/s  for  InAs.  Because  the  0.02  ML/s  growth  rate  for  GaAs  growth  is  too 
high  to  distinguish  the  dots. 

7)  High  resolution  transmission  electron  microscopy  study  of  InAs  on  Si  (OOP  : 

The  detailed  structure  of  InAs  on  Si(001)  is  studied  via  high  resolution  TEM.  The  strain 
is  relieved  by  introducing  misfit  dislocation  at  very  early  growth  stages,  with  island  size 
larger  than  4nm,  as  shown  in  Fig.  31.  Both  glissile  60  0  dislocation  and  sessile  90  0 
dislocations  are  observed  by  TEM.  At  the  beginning  of  dislocation  formation,  the  60  0 
dislocation  is  formed.  With  increasing  size  of  island,  the  90°  dislocations  appear.  It  is 
very  likely  that  the  90°  dislocations  are  formed  by  two  60°  dislocation  with  increasing 
dislocation  number,  as  shown  in  Fig.  32 


Figure  31.  60°  misfit  dislocations  formed  at  the  interface  (indicated  by  red  lines). 


27 


<no> 


hA  *11  i. 

Ja 

>!v^  Ay 


Figure  32.  Both  60°  and  90°  dislocations  formed  at  the  interface  between  Si  and  InAs.  The  enlarged  images 
show  the  detailed  atomic  structure  of  dislocations. 

Summary 

We  have  demonstrated  the  key  parts  of  a  Si-based  laser  structure  with  optimized 
design.  Relatively  low  propagation  loss  was  achieved  even  in  a  SiGe  waveguide  with 
metal  cladding,  high  injection  current  of  more  than  2  kA/cm2  can  be  easily  achieved  in 
such  a  low  loss  waveguide,  and  single  mode  resonance  was  achieved  in  a  SiGe  DBR 
resonant  cavity.  In  addition,  we  have  carried  out  theoretical  calculation  for  coupled 
SAQD  systems  and  formulated  the  distinct  advantage  of  a  coupled  SAQD  system  for 
laser  applications.  We  also  achieve  controllable  growth  of  InAs  quantum  dots  on  silicon, 
and  a  detail  study  of  structure  characteristics  of  such  quantum  dots  was  carried  out. 

The  crucial  finding  from  this  program  is  that  the  very  small  critical  size  of  InAs 
SAQD  for  dislocation  and  the  associated  extreme  quantum  confinement  show  that  there 
might  not  be  a  window  in  the  dot  size  for  which  electron-hole  pairs  can  be  confined  and 
the  dots  are  coherently  strained.  It  can  be  used  as  a  guide  for  future  funding  decisions  on 
the  subject  of  lasers  on  Si. 
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